Astronomy & Astrophysics manuscript no. 


February 2, 2008 


(DOI: will be inserted by hand later) 





(N 
O 
O 
(N 

< 

O 

(N 



> 
(N 

m 
oo 
o 

(N 
O 

0,; 

6 

c3 



Quasi-molecular lines in Lyman wings of cool DA white dwarfs 
Application to FUSE observations of G 231-40 

G. Hebrard 1 , N. F. Allard 1 ' 2 , I. Hubeny 3 , S. Lacour 4 * 1 , R. Ferlet 1 , and A. Vidal-Madjar 1 



1 Institut d'Astrophysique de Paris, CNRS, 



boulevard Arago, F-75014 Paris, France 



2 Observatoire de Paris-Meudon, LERMA, F-92195 Meudon Principal Cedex, France 

3 NOAO, 950 North Cherry Avenue, Tucson, AZ 85726, USA 

4 Department of Physics and Astronomy, Johns Hopkins University, Baltimore, MD 21218, USA 
Received accepted ... 

Abstract. We present new theoretical calculations of the total line profiles of Lyman a and Lyman f3 which 
include perturbations by both neutral hydrogen and protons and all possible quasi-molecular states of H2 and 
Hj. They are used to improve theoretical modeling of synthetic spectra for cool DA white dwarfs. We compare 
them with FUSE observation of G 231-40. 

The appearance of the line wings between Lyman a and Lyman j3 is shown to be sensitive to the relative abundance 
of hydrogen ions and neutral atoms, and thereby to provide a temperature diagnostic for stellar atmospheres and 
laboratory plasmas. 
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1. Introduction 

Structures in the Lyman a and Lyman (3 line wings 
have been identified with free-free transitions which take 
place during binary close collisions of the radiating H 
atom and a perturbing atom or ion (Allard et al. 1998a, 
1998b, 1999). The characteristics of these features (posi- 
tion, amplitude, and shape), due to the formation of quasi- 
molecules during collisions between the radiating atom 
and perturbers, depend directly on the potential energy 
curves correlated to the atomic levels of the transition 
(Allard & Kielkopf 1982). 

Two satellite absorption features at 1058 A and 
1076 A due to collisions of atomic hydrogen with pro- 
tons were first identified in the spectrum of the DA white 
dwarf Wolf 1346, as observed with the Hopkins Ultraviolet 
Telescope (Koester et al. 1996). These satellites in the red 
wing of Lyman j3 are in the Far Ultraviolet Spectroscopic 
Explorer (FUSE) spectral range (Moos et al. 2000); fur- 
thermore, Lyman f3 profiles are also the subject of an on- 
going study of the far ultraviolet spectrum of dense hy- 
drogen plasmas. 

In Allard et al. (1998a) we presented theoretical pro- 
files of Lyman /3 perturbed solely by protons. The calcula- 
tions were based on the accurate theoretical molecular 
potentials of Madsen & Peek (1971) to describe the inter- 
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action between radiator and perturber, and dipole transi- 
tion moments of Ramaker & Peek (1972). The line profiles 
were included as a source of opacity in model atmospheres 
for hot white dwarfs, and the predicted spectra com- 
pared well with the observed ORFE "US and FUSE spectra 
(Koester et al. 1998; Wolff et al. 2001). 

Ab initio calculations of Drira (1999) of electronic tran- 
sition moments for excited states of the H2 molecule and 
molecular potentials of Detmer et al. (1998) allowed us to 
compute Lyman (3 profiles perturbed by neutral atomic 
hydrogen (Allard et al. 2000). The appearance of a broad 
satellite situated at 1150 A makes necessary to take into 
account the total contribution of both the Lyman a and 
Lyman (3 wings of H perturbed simultaneously by neutrals 
and protons. 

We show that the shape of the wings in the region 
between Lyman (3 and Lyman a is particularly sensitive 
to the relative abundance of the neutral and ion perturbers 
responsible for the broadening of the lines. 

These new profiles have been used to predict synthetic 
spectra for cool DA white dwarfs which present structures 
at 1600 A and 1400 A in the Lyman a wing due respec- 
tively to quasi-molecular absorption of the H2 and 
molecules. These last two structures have been demon- 
strated to be a very sensitive temperature indicators in DA 
white dwarfs. The relative strength of these two satellite 
features depends very strongly on the degree of ionization 
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in the stellar atmosphere, and thus on the stellar param- 
eters T e g and logg (Koester & Allard 1993; Koester et 
al. 1994; Bergeron et al. 1995). 

2. Theoretical line profiles 

2.1. Theory 

We use a general unified theory in which the electric dipole 
moment varies during a collision; a detailed description of 
the theory as applied to the shape of the Lyman lines 
has been given by Allard et al. (1999). The obtained line 
profiles fit the spectra of laser-produced hydrogen plasmas 
(Kielkopf& Allard 1998). 

Our approach requires prior knowledge of accurate the- 
oretical molecular potentials to describe the interaction 
between radiator and perturber, and knowledge of the 
variation of the radiative dipole moment with atom-atom 
and atom-ion separation for each molecular state. This ef- 
fect is important when the dipole moment varies in the re- 
gion of inter-nuclear distance where the satellite is formed, 
and thus cannot be neglected. In the case of Lyman (3 
satellites, due to H-H + collisions, we have shown that 
large changes (up to 60 %) in the intensity of the satel- 
lites may occur when the variation of the dipole moment is 
taken into account (Allard et al. 1998a). This result is also 
valid for other lines; it increases by a factor of about 2 the 
main satellites of Lyman a (Allard et al. 1999). Previous 
line profile calculations using constant dipole moment 
have been used to interpret IUE (International Ultraviolet 
Explorer) and Hubble Space Telescope spectra by Koester 
& Allard (1993), Koester et al. (1994), and Bergeron et 
al. (1995). The synthetic spectra presented here used im- 
proved Lyman a line profiles of Allard et al. (1998b) which 
have been already included in stellar atmosphere programs 
for the computation of stellar atmosphere model and syn- 
thetic spectra of A Bootis stars. A comparison of these 
calculations with observations made with the IUE demon- 
strated that these last improvements are of fundamental 
importance for obtaining a better quantitative interpreta- 
tion of the spectra and for determining stellar atmospheric 
parameters (Allard et al. 1998b). 

2.2. Lyman (5 in H-H collisions 

The Lyman profiles and satellites are calculated at the 
low densities met in the atmospheres of stars. The typical 
particle densities (10 15 to 10 17 cm -3 ) allows us to use an 
expansion of the autocorrelation function in powers of den- 
sity as described in Allard et al. (1994) and Royer (1971). 
Line profiles are normalized so that over 1/A = u> (cm -1 ) 
they integrate to 1. 

The only line feature of the Lyman (3 profile is a broad 
absorption line satellite situated at 1150 A due to the 
B"B X S+-X X E+ molecular transition of H 2 (Fig.|). The 
recent ab initio calculations of A. Spielfiedel (2001, private 
communication) have shown that for the isolated radiat- 
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Fig. 1. Lyman (3 satellite due to H-H collisions. The neu- 
tral density uh is 1 x 10 16 cm~ 3 . 

ing atom ( R — > oo ) this transition is not asymptotically 
forbidden as it was explicitly stated in Allard et al. (2000). 

The line satellite shown in Fig. [l] presents a shoul- 
der at 1120 A, a similar shape has been obtained for 
the 1600 A satellite. In Fig. 6 of Allard et al. (1999) 
both theory and experiment show an oscillatory struc- 
ture between the satellite and the line, with a minimum 
at about 1525 A. These oscillations are an interference 
effect (Royer 1971; Sando & Wormhoudt 1973), and are 
expected to depend on the relative velocity of the collision 
and therefore on temperature. 

This satellite of Lyman (3 is quite far from the un- 
perturbed Lyman (3 line center, actually closer to the 
Lyman a line. It is therefore necessary to take into ac- 
count the total contribution of both the Lyman a and 
Lyman (3 wings of H perturbed simultaneously by neu- 
trals and protons and to study the variation of this part 
of the Lyman series with the relative density of ionized 
and neutral atoms. 

2.3. Lyman (3 in H-H and H-H+ collisions 

In Allard et al. (1998a) we presented Lyman (3 profiles 
perturbed by protons. The line profile calculations were 
done without using the expansion in density and then were 
valid from the center to the far wing and allowed a com- 
parison of the amplitudes of the satellites to the line core. 
The profiles at different densities of H + were displayed 
in Fig. 8. 

In Fig. ^ we show the sum of the profiles of Lyman a 
and Lyman (3 perturbed by collisions with neutral hydro- 
gen and protons for different densities nu ■ 

We can see that a ratio of 5 between the neutral and 
proton density is enough to make the quasi-molecular H 2 
satellite appear in the far wing. Note that the 1120 A 
shoulder of the satellite is still visible in the total profile. 
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Fig. 2. Total profile of Lyman a and Lyman (3 perturbed 
by neutral hydrogen and protons. Three different neutral 
densities (1 x 10 17 , 5 x 10 16 , 1 x 10 16 cm' 3 , top to the 
bottom) are compared for a fixed ion density (1 x 10 16 
cm~ 3 ). 



3. Synthetic spectra for DA white dwarfs 

We now take into account these absorption features of the 
H+ and H 2 quasi-molecules to calculate synthetic spectra 
for DA white dwarfs. Atmosphere models and synthetic 
spectra have been calculated using the computer programs 
TLUSTY and SYNSPEC (Hubeny 1988; Hubeny & Lanz 
1992, 1995). The LTE model atmospheres assume a pure 
hydrogen composition including the quasi-molecular opac- 
ities. 

Below T c ff=15 000 K, the atmosphere of DA white 
dwarfs become convective, and model atmospheres for 
these stars are usually calculated using several variants of 
the standard mixing-length theory. According Bergeron et 
al. (1992), we used the so-called ML2 prescription with 
a = 0.6, which was shown to provide excellent internal 
consistency between ultraviolet and optical temperatures. 

The temperature range where both H2 and Hj 
Lyman (3 satellites are visible is relatively small; roughly 
between 14500 to 11000 K (Fig. |-upper panel). This de- 
pendence of the theoretical profiles on temperature is ex- 
tremely strong because of the relative importance of per- 
turbations by neutral versus ionized hydrogen. The satel- 
lite appearance is then very sensitive to the degree of ion- 
ization and may be used as a temperature diagnostic. 

In order to see all the Lyman line satellites, synthetic 
spectra are plotted from 1020 to 1800 A (Fig. ||dower 
panel). Below T c s= 13000 K the H^ starts to disappear 
but the broad H2 satellite at 1150 A is still present and 
the one at 1600 A slowly increases. 

4. FUSE observations of G 231-40 



4.1. Observations and data reduction 

G 231-40 (WD 2117 + 539) is a cool DA white dwarf lo- 
cated at Galactic coordinates I = 95.0° and b = —3.3°; it is 
a relatively bright white dwarf: V — 12.3. This target was 
observed with FUSE as part of the French Guaranteed 
Time Observing Programs (program Q210). Two expo- 
sures were obtained on 2001 July 3 in time-tagged pho- 
ton address mode (TTAG) with the object in the large 
aperture (LWRS). The total duration was ~ 1.6 h (see 
Table 0). Details of the FUSE instrument may be found 
in Moos et al. (2000) and Sahnow et al. (2000). 

The one-dimensional spectra were extracted from the 
two-dimensional detector images and calibrated using ver- 
sion 2.0.5 of the CalFUSE pipeline. The eight FUSE de- 
tector segments of the two exposures (i.e. 16 spectra) were 
co-added and projected on a 0.16 A-pixel base, i.e. pixels 
about 25 times larger than the original FUSE detectors 
pixels. This degradation of the FUSE spectral resolution 
(typically A/AA ~ 15000 for this kind of target in the 
large slit; see Hebrard et al. 2002; Wood et al. 2002) is of 
no effect on the shapes of the large stellar features which 
we study and allows us to increase the signal-to-noise ra- 
tio. At this resolution, no spectral shifts were detected 
between the 16 different co-added spectra. 

Poor quality edge of each segment were not included 
in the sum. We used the following spectral ranges for each 
segment: 1010 - 1085A (SiCIA), 910 - 992.5 A (SiClB), 
922 - 1005 A (SiC2A), 1016 - 1101 A (SiC2B), 1004 - 
1081.5 A (LiFlA), 1096 - 1187 A (LiFIB), 1087 - 1181 A 
(LiF2A), and 1012 - 1073 A (LiF2B). As the LiFIB seg- 
ment presents a known large-scale distortion ( "the worm" , 
D. Sahnow 2000, private communication) in the flux cal- 
ibration around 1120 — 1170 A (see, e.g., Hebrard et 
al. 2002), we normalized the large-scale shape of the LiFIB 
segment to the shape of the LiF2A segment, using a poly- 
nomial fit. We also did not include the part of the LiF2A 
spectra in the ranges 1134.6—1135.6 A and 1151.5—1153 A 
because they were deteriorated by the so-called "walk" 
problem, which is a distortion in the FUSE spectra at air- 
glow wavelengths (A 1134 A Ni and A 1152 A Oi lines 
in the present case) caused by variation of position with 
pulse height (D. Sahnow 2001, private communication). 
The final FUSE spectrum is plotted in Fig. |. Note that 
the interstellar Ni A 1134 A triplet is detected on this line 
of sight. 

Since 2001, repeated observations of standard white 
dwarfs indicate a slow degradation in the effective area of 
the FUSE spectrograph. The decline affects the flux up 
to 20%, and that mainly for the LiF2A channel. A dime- 
dependent flux correction is in development for a future 
version of the pipeline, but is not included in the CalFUSE 
version 2.0.5. 

In order to complement the 910 — 1187 A FUSE spec- 
trum, we retrieved two IUE low dispersion spectra of 
G 231-40 from the STScI MAST archive. These spectra 
were obtained in March 1983 and they were reduced us- 
ing the NEWSIPS processing (see Table |). They yield the 
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Fig. 3. Synthetic spectra for DA white dwarfs including Lyman a and Lyman [3 satellites, on spectral ranges 1020 - 
1200 A (upper panel) and 1020 - 1800 A (lower panel). T c ff=14500, 14000, 13 500, 13 000 K (from top to bottom) 
and log g = 8.0. The flux, at the stellar surface, is integrated over angles. 



Table 1. Observation log 



Instrument 


Observation 


Date of observation 


Aperture 


Exposure time (sec) 


Spectral extraction 


FUSE 


Q2100101001 


2001 July 3 


LWRS (30" x 30") 


2940 


CalFUSE 2.0.5 


FUSE 


Q2100101002 


2001 July 3 


LWRS (30" x 30") 


2864 


CalFUSE 2.0.5 


IUE 


SWP19532 


1983 March 24 


SWLA (10" x 20") 


720 


NEWSIPS 


IUE 


SWP 19533 


1983 March 24 


SWLA (10" x 20") 


720 


NEWSIPS 



spectral coverage 1150 - 1980 A at a resolution of about 
7 A (Wegner 1984). The flux calibration of the FUSE 
and the IUE spectra are in good agreement (see Fig. |], 
lower panel). 

4.2. Comparison with theoretical spectra 

As a change in temperature can be compensated by 
a change in surface gravity we have used logg deter- 
mined from optical analysis. Holberg et al. (1998) give 
(14490; 7.85) for (T c g; logg), the atmospheric parame- 
ters being determined from spectroscopic observations by 
Bergeron et al. (1992). Model atmospheres and synthetic 



spectra were calculated in the range 15 500 to 14 500 K by 
steps of 100K, using log<7=7.85. 

H opacity (Wishart 1979) has been added in the at- 
mosphere model calculation and gives a feature at 1130 A 
also present in the FUSE observation (see Fig. ^). 

Fig. ^- upper panel shows the comparison of the FUSE 
observation with our best fit obtained for T e ff=14800 K. 
The decre asing of sensitivity on the segment LiF2A 
(Sect. 4T) is the main explanation of the difference be- 
tween the data and the model between 1120 and 1180 A. 
The H 2 satellite (1150 A) is no more visible, as expected, 
in this range of temperature. 



G. Hebrard et al.: Quasi-molecular lines in Lyman wings of cool DA white dwarfs 



5 



6 

4 h 

2 



^ 

CD 




1050 



1 00 



1 1 50 



1 200 




200 



1 400 
■'avelength 



600 



1 800 



in A 



Fig. 4. FUSE (1020 - 1187 A, upper and lower panels) and WE (1150 - 1800 A, lower panel) spectra of G 231-40 
compared with a theoretical model for T c ff= 14800 K and logg = 7.85. Emissions at Lyman a, Lyman (3, and 1027.5 A 
are due to Hi and Oi airglow. The difference between the data and the model between 1120 and 1180 A is mainly 
due to an uncorrected decreasing of sensitivity on the FUSE segment L1F2A (see Sect. 4.1). Notice the H~ feature at 
1130 A. 



Fig. |] shows a good agreement between the FUSE and 
WE observations of Lyman lines and predicted spectra 
over the whole range from 1020 to 1800 A using atmo- 
spheric parameters determined from the optical range. 

The feature near A 995 A due to a Lyman 7 satellite, 
which is visible in the spectra of hotter objects, is no 
more present in the temperature range of G 231-40 (the 
flux of G 231-40 is below 1 x 10" 14 erg cm~ 2 s" 1 A" 1 
at 995 A). It was visible in FUSE observations of 
CD -38° 10980 (Wolff et al. 2001), HUT spectrum of 
Wolf 1346 (Koester et al. 1996), and in some ORFEUS 
spectra (Koester et al. 1998). It was also recently seen in 
FUSE spectra of Sirius B (Holberg et al. 2002). This fea- 
ture, due to a Lyman 7 satellite, is now included in our 
calculations. We have shown that a larger temperature is 
required to get it detectable (Allard et al. 2002). 



5. Conclusion 

From theoretical profiles including the new Lyman [3 opac- 
ities it can be predicted that Lyman (3 satellites should be 
detectable roughly between 25000 K to 11000 K. Within 
this range in a very small domain of temperatures, the 
whole profile should present H2 and Hj satellites of the 
Lyman a and Lyman (3 lines. FUSE observations of cooler 
white dwarfs in the ZZ Ceti range stars would offer the 
best opportunity to determine accurate stellar parameters 
as T e g and log g for these stars. 
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